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Abstract 

For a particular family of long-range potentials V, we prove that the eigenvalues of the 
indefinite Sturm-Liouville operator A = sign(a:)(—A + V(x)) accumulate to zero asymptotically 
along specific curves in the complex plane. Additionally, we relate the asymptotics of complex 
eigenvalues to the two-term asymptotics of the eigenvalues of associated self-adjoint operators. 


1 Introduction 


Given a real-valued potential V such that 


V € L°°(R), lim V[x) = 0, 

rr—>-d=oo 


lint sup x 2 V(x) < - 
X —>-±00 4 


( 1 ) 


consider a one-dimensional Schrodinger operator in L 2 (M) 

T:=TV: =-;S +F W’ 

Dorn(T) := {/ € L 2 {R) \ fj'e AC(R), Tf € L 2 (M)} . 


( 2 ) 


It is well known that in this case the spectrum Spec(T) is bounded from below, the essential 
spectrum Spec ess {T) = [0, oo), and the negative spectrum Spec(T)n(—oo, 0) consists of eigenvalues 
accumulating to zero from below. 

Let J := sign(x) be the multiplication operator by ±1 on M±. In what follows we consider the 
point spectrum of the operator 


A := A v := JT V , Dom(A) = Dom(T). (3) 

This operator is not self-adjoint (and not even symmetric) on L 2 (R), and its spectrum need not 
therefore be real. However, as J* = J= J, A can be treated as a self-adjoint operator in the 
Krein space (L 2 (M), [•, ■]) with indefinite inner product 

[f,g\ ■= (Jf,g) L ( R) = J f(x)g(x) sign(x)dx 
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or equivalently as a J-self-adjoint operator Azlo . Operators of type ([ 3 ]) have been studied both 
in the framework of operator pencils, cf. |DaLel [Maj. and of indefinite Sturm-Liouville problems 
|BeKaTrl [BiTrl [RaTrl [La]. 

In both settings the literature is extensive, starting mostly with Soviet contributions in the 
1960s, including those by Krein, Langer, Gohberg, Pontryagin and Shkalikov. We refer to {MHO] 
for reviews and bibliographies. In particular due to its many applications, for example in control 
theory, mathematical physics and mechanics, the field is still very active, with recent works on 
the theoretical,as well as numerical, aspects, (see e.g. DaLe EILePol IHiTrVDl l|Vej and references 
therein). 

In the special case of indefinite Sturm-Liouville operators, it is well known that for positive 
potentials, V > 0, the spectrum of Ay is real and the operator Ay is similar to a self-adjoint 
operator CuLa. CuNa. Py| IKoj . At a very basic level, this can be seen from the following abstract 
construction: if R and S are self-adjoint operators with R > 0, then, under mild restrictions, the 
spectrum of R~ l S is the same as the spectrum of the self-adjoint operator R^ 1 / 2 SR~ 1 / 2 , and is 
therefore real. 

The case V G L 1 (M, (1 + |x|)dx) has been considered in {KaKoMaf . where it is shown that Ay 
is self-adjoint iff Spec(Ay) = M. Finally, for (quasi-)periodic finite zone potentials, KaMa explores 
some conditions under which Ay is similar to a self-adjoint or a normal operator. For a review of 
indefinite weighted Sturm-Liouville problems on a finite interval, see fFTj. 

Let us return to our original problem ([ 3 ]). Recently, there was a rapid growth of interest in 
the case of non-positive potentials, especially by Behrndt, Trunk, and collaborators [ Be07, [BeKaTrl 
BeTrl IBePhTrj . clarifying the structure of their spectra and other properties as well as stating new 
conjectures on rather unusual spectral behaviour [ Bel3 j. 

The following known results are a particular case of [ BeKaTr , Theorem 1 and Theorem 2] and 
[BePhTrl Theorem 4.2]. 

Proposition 1 . For the operator A in ([ 3 ]), 

(a) Spec(A) is symmetric with respect to M. 

(b) Spec ess (A) = R. 

(c) Spec(A) \ M consists of eigenvalues of finite multiplicity. 

(d) No point ofM \ {0} is an accumulation point of non-real eigenvalues of A. 

(e) At least one of the following statements is true: 

(i) The non-real eigenvalues of A accumulate only to 0; 

(ii) There exist embedded eigenvalues of A in M + that accumulate to 0; 

(iii) There exist embedded eigenvalues of A in M_ that accumulate to 0; 

(iv) The growth of A i-)- {A — p)~ 1 near zero is not of finite order. 

(0 If additionally V is even, V(x) = V(—x), then Spec(A) is also symmetric with respect to iM. 


(g) The non-real spectrum of A is contained in the strip | Irn p\ < 21| ||oo. 

Despite the amount of information on the structure of the spectrum of A, known proofs of 
Proposition |T] are not constructive and, in fact, we do not even know a priori which of the four 
statements (e)(i)-(iv) are true for a particular given potential. Some numerical experiments, cf. 
[ BeKaTr ], have recently led to conjecture that statement (e)(i) in Proposition 1 may hold for many 
potentials satisfying Q, see Figure[l] 

In this paper we prove that for a particular family of potentials 

V(x) = V 1 (x) = - T ^— 7 > 0, (4) 
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Figure 1: A numerical example showing accumulation to 0 of complex eigenvalues (red diamonds) 
of the operator A 7 , 7 = 2.5. The magenta and white circles on the negative real axis are the 
eigenvalues of T 7 corresponding to the eigenfunctions which are even or odd with respect to zero, 
cf. |BeKaTr| . 

Proposition l(e)(i) holds. Moreover we also prove (Theorem [ 6 ]) that the complex eigenvalues of 

T = T 7 := T Vi 

accumulate to zero asymptotically along specific curves in the complex plane, and that the explicit 
asymptotics of complex eigenvalues of T 7 can be obtained from the asymptotics of eigenvalues of 
the self-adjoint operator 

A = A 7 := Ay 1 (5) 

(or, more precisely, from the eigenvalues of its restriction on either even or odd (with respect to 
zero) subspace). We also extend these results to the more general non-symmetric potentials 

f-—^ if x > 0 

V 1 -, 1+ {x) = l 7 ]^ 7+>7-f®+' (6) 

- , , , if x < 0 

l l + \x\ 

The rest of the paper is organised as follows. Section 2 contains the statements of our main 
results. The proofs, as well as some numerical examples, are in Sections 3-5; they are based on the 
explicit expressions for Jost solutions of the differential equation 

d 2 7 

on M + and a rather delicate asymptotic analysis involving Kummer functions. A brief exposition 
of some auxiliary results, mainly due to Temme [TeJ, which we use in our proofs, is given in the 
Appendix. 

2 Sharp asymptotics of the eigenvalues of the self-ajdiont operator 

Let T and denote the restrictions of the operator T 7 to M+ with Dirichlet and Neumann 
boundary condition at zero, resp. By the spectral theorem, for symmetric potentials V 7 (x) 

Spec(T 7 ) = Spec (T^) U Spec(T^) 
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with account of multiplicities. Let — Ajf (7) denote the eigenvalues of T#, # = D or N, ordered 
increasingly. In what follows we often drop the explicit dependence on 7. 

It is well-known that — A# < 0 and —A^ < —A^ < —A^ +1 for all n G N, and also that 

—A^f -7 0— as n -7 00 . 

Before stating our main results, we need some additional notation. 


Definition 2. Let F denote the class of piecewise smooth functions F : M_|_ -7 M which have a 
discrete set of singularities (with no finite accumulation points). At each singularity both one-sided 
limits of F are ±00 and differ by sign. Assume for simplicity that 0 is not a singularity of F, and that 
-F(O) = 0. For F £ F we denote by Qp(x) the continuous branch of the multi-valued Arctan(i ? (x)) 
such that ©^(0) = 0. 

Remark 3. Away from the singularities of F, the function @f(x) can be written in terms of the 
ordinary arctan(F(x)) (which takes the values in [—§,§]) and the total signed index of F on [0, x], 
which we denote by Zp(x), and which is defined as the total number of jumps from +00 to —00 
on [0,x] minus the total number of jumps in the opposite direction: 


Zf(x ) 


E 


E 


{rG(0,fc]| lim F(t)=-\-oo} {t(E(0,:c]| lim F(t)=— 00 } 

£—ft — t —ft — , 


(7) 


Then 

&f(x) = arctan (F(x)) + itZf(x). (8) 

Obviously, Q_f(x) = —@f(x). 

Our first result gives sharp two-term asymptotics of eigenvalues (accumulating to zero) of the 
self-adjoint operators Tf. 

Theorem 4. As n -7 oo, 


A ?W-^( 1 -S e "‘W + 0 (?))- 

A ” W = £( 1 -7 e *»W + 0 (^))’ 


where 




Jk (2^7) 
Y k ( 2Vt)’ 


and Jk and Yj. denote the Bessel functions of the first and second kind, respectively. 
This immediately implies 


Corollary 5. As n —> oo, 

A£ - A % = 0 ((A#) 3 / 2 ) . 

Some old papers found both in physical and mathematical literature addressed the problem of 
approximating the eigenvalues of Schrodinger operators with shifted Coulomb potentials, see e.g. 
Gi 0 EQ and references therein. However they were considering somewhat different asymptotic 
limits, and to the best of our knowledge the two-terms asymptotics of Theorem [4] are new. 
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Figure 2: Rii'y) and i?o( 7 ) (resp. magenta and blue line) from Theorem [4] and the corresponding 
@^( 7 ) and Qr 0 ( 7 ) (resp. magenta and blue dashed line). 



Figure 3: Approximate eigenvalues of T 7 , 7 = 2.5, as described in Theorem [4j The magenta (resp. 
white) circles correspond to —A® (resp. —A^). The solid (resp. dashed) lines are the left-hand 


sides of (25) (resp. (26)) whose roots are the eigenvalues. Only curves and eigenvalues up to —0.01 


are displayed. 

3 Sharp asymptotics of the eigenvalues of the non-self-ajdiont op¬ 
erator 

Our main result is the following 

Theorem 6. (i) The eigenvalues of A 7 lie asymptotically on the curves 

|Im p\ = T(7) I Rep \ 3 / 2 + O ((Re /u) 2 ) , p 0 , 

where 


and 


T(7) = — log ^ , 
7T7 1 + <r( 7) 


q( 7) : = M 2 Vt) + y o{ 2 Vt) Y i( 2 Vr)) ■ 


(9) 

( 10 ) 

(ii) 


(ii) More precisely, the eigenvalues {p} n & n of in the first quadrant (ordered by decreasing real 
part) are related to the absolute values A# of the eigenvalues of the self-adjoint operators T#, 
# = D or N, by 

p n = A° + T-(7)(Af ) 3 / 2 + O ((A^) 2 ) = A^ + T + ( 7 )(A^) 3 / 2 + O ((A^) 2 ) (12) 

as n —> 00, where 

TT( 7 )=A arctm (_J_). (is) 

The expressions for eigenvalues in the other quadrants are obtained by symmetries with respect 
to M and iM. 
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Before proceeding to the proofs, we want to discuss the statements of Theorem [6] in more detail. 


Remark 7. (a) It is immediately seen from (10) and (13) that 


T( 7 ) = Im T _ ( 7 ) = Im X + ( 7 ). 


(b) If we introduce two functions r T : M + — > C by 

T^(t) = t + r^t 3 / 2 , 


then 

IniT~(f) = IniT + (t) = T f 3//2 , 

which is just another way of writing (a). (We have dropped the dependence on 7 here for 
clarity.) 


(c) The statement (ii) of Theoremj6l contains in fact several results. Taking the imaginary parts 
of (12) leads immediately to ([£]) by Corollary [5] The equalities obtained by taking the real 
parts (12) are more intricate: they indicate that, up to the terms of order (Ajf) 2 , the values 
of Re(r _ (A^)) and Re(T + (A^)) coincide. 


In other words, we can construct the eigenvalues of the non-self-adjoint operator A 7 by per¬ 
turbing either the anti-symmetric self-adjoint eigenvalues A^ or the symmetric self-adjoint 
eigenvalues A^, and the different formulae still lead to the same result, modulo higher-order 
terms. 


(d) The only previously known bound, see Proposition 1(g), implies only that for our potential V 1: 

| Im /i| < 2 7 . 

The typical eigenvalue behaviour is illustrated in Figure[4] 



Figure 4: The red diamonds are numerically computed exact eigenvalues of Ay, 7 = 2.5, lying in the 
first quadrant. The complex parametric curves fi = r~(t) (the magenta solid line) and /i = r + (f) 
(the blue dashed line) are as in Remark [TJb). The approximated complex eigenvalues, computed by 
the first part of formula (12) are shown as magenta squares, and computed by the second part of 
formula (12) are shown as white squares. The absolute values of the eigenvalues of T 7 are marked on 
the real line in the same way as in Figure [3] The dotted arrows are to indicate which real eigenvalue 
“produces” the corresponding complex ones. 
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4 Explicit form of the Jost solutions 


4.1 Solutions of the equation on the half-line 

Let /c be a general solution of the differential equation 

-j-2f(y)--f(y) = £f(y)> y^ R +- 

ay y 

The change of variables 

V = + 

relates /g with solutions of the differential equation 

d 2 7 

~ l6R +’ 


by 


9n,d x ) = fvh 2 ( 7 (® + 1 )) • 


(14) 

(15) 

(16) 

(17) 


4.2 Explicit solutions of the differential equation (16) 

Assuming £ E C+, we will write 

s := -V=£ (18) 

where we take the principal branch of i.e. the uniquely determined analytic branch that maps 


M_l_ into itself. Obviously s 2 = £. 
With the ansatz 

equation (p~4]) can be reduced to 


f(y) = ye lsy h{y ), 


d 2 d 

~ V dyi h ^ ~ ~ 2[s y^ h ^ + ( 2is ~ !) %) = 0. 


(19) 


( 20 ) 


With the additional change of variables 

w = 2i sy (21) 

we arrive at a particular case of the Kummer Hypergeometric Equation [ OlLoBoC I, Chapter 13.1] 


d 2 ~ d ~ 

w- — ~h(w ) + (b — w)- — h(w) — ah(w) = 0 
aw z aw 


( 22 ) 


with 


° = 1 _ 2ls’ 


= 2 , 


w = 2i sy, h{w) = h f 


The two linearly independent solutions of (22) are known as Kummer hypergeometric functions 

M(a,b;w) and U(a,b;w) so that 

h(w) = CiU(a , b ; w) + C 2 M(a, 6; w), C\.C 2 = const. 

Hence the solutions f^(y) of ( p~4] ) are of the form 

fdv) = ye~ isy ( CiU (l - ^,2; 2i sy^j + C 2 M (l - ^,2; 2isy)) . (23) 

Thus, by ([TT]), the solutions of ( fl6| ) are of the form 

Y (x) = 7(1 + x)e~^ 1(1+x) (ciU ^1 - 2 ^— , 2 ; (1 + *)) 


9A,7V 


+ C 2 M 1 - 


2\/—A 


, 2; 2 \J~n (1 + x)^ ^ . 
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4.3 The Jost solutions of ( |16| ) 

It is well known, see e.g. ( OlLoBoCI . 13.7.1 and 13.7.2], that the first order asymptotic behaviour 
of the Kummer Hypergeometric Functions is given, as |w| —y oo, by 


U(a, 6; w) ~ w 
M(a, b ; w 


e w w b a 


+ 


r(a) T(6 — a) ’ 


37t 37r 

~Y <ar s( w ) < y> 

7 T 3 tT 

~2 < ai 'g(w) < y, a,b-a 0 -N U {0}, 


where T(-) stands for the usual Gamma function. 

For fj, G C\M + , we have — i/yZ C {z G C | Re z < 0}, and therefore 


£7(1 — 


7 


2yYY 


2 ; 2y^PY) ~ (2 


-1 


and 


Ml- 


^=,2; W=S) ~ 


as y — y oo. 

This in turn implies that the /c and gy defined above are L 2 (M + ) if and only if C 2 = 0. For 
convenience, we choose further on the normalisation C\ = 1. 

The L 2 (M_|_) solutions of (16) are called the Jost solutions. We denote them by 

:= ^,7(^)101=1,02=0 = 7(1 +x)e'^^ +x) U (l- 2 ^— ,2; 2(1 + x )y/=jij . (24) 

5 Proof of Theorem 0] 

5.1 The characteristic equations for the self-adjoint problem 

It is well known that T has a negative discrete spectrum accumulating to zero. We need a more 
detailed knowledge of the asymptotics of its eigenvalues and of the corresponding eigenfunctions. 

It follows from the arguments of the previous sections that up to a scaling constant the eigen¬ 
functions of the self-adjoint problem (16) are given, on M+, by 

V>a(7, x) = 7(1 + x)e _v/X(1+ - T) U 1 - -^=,2] 2(1 + s)\/Aj . 

The eigenvalues —X^ of the self-adjoint operator with Dirichlet boundary conditions at zero 
are thus given by the solutions of i/j\( 7 ,0) = 0, i.e. 


7e 


-7a 


2V 7 ! 


U - 


7 


2^X 


0,2VA ) = 0 


(25) 


The eigenvalues — X^ of the self-adjoint operator with Neumann boundary conditions at 
zero are given by the solutions of 4zip\{j,x) 


x=0 


7 

= 0, i.e. 


7e 


-7a 


A 2 


((7 - 2a/A)£7 (-^,-1,2^) + 2a/A(TA - 1)U (-^=, 0, 2>/a)) = 0. 


(26) 


The solutions of transcendental equations (25) and (26) can be computed numerically, although 
it is a non-trivial task as the left-hand sides of these equations oscillate wildly for small A (cf. Figure 
[6]). We instead use asymptotic techniques to approximate the Kummer functions as A —>• 0 and to 
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control their oscillations. A quick took at (25) and (26) shows that we require asymptotic formulas, 
as A —> 0+, for 


U - 


(-^,c ; 2Va), c £ {0, —1}, 


(27) 


Unfortunately, it is a difficult task — the corresponding formulas, are not, in fact, in the standard 
references. We rely, instead, on the results from the forthcoming book [Tej which we summarise 
and adapt in the Appendix. 

5.2 Asymptotic solutions of a transcendental equation 

A crucial element of our analysis is the investigation of the large /-c-roots of the equation 

tan(7«) = G(k, 7) ( 28 ) 

where 7 is treated as a parameter, and where G depends analytically on k in the vicinity of k = 00 
and, to leading order, is of class T as a function of 7. The required results are summarised in the 
following 

Lemma 8. Let G(k, 7) be an analytic function of k around k = 00 such that 
G(rt,7) = G 0 (7) (l + 0 (k -1 )) , as k- t 00, 

Gq £ fF, and the O terms are regular in 7. Then the solutions n n ( 7), ordered increasingly, of the 


equation ( 28 ), are given, as n 00 , by 


«n(7) = — + -Qgo(t) + °( n X )- 
7 7 


(29) 


The proof of Lemma [8] is in fact immediate as soon as we recall Definition [2] of © and the fact 
that tan is 7r-periodic. 

Considering additional terms in the expansion of G one can get additional terms in the expansion 


of k„. This is in fact what we do in more detail in Section 6.2 


5.3 Approximation of Dirichlet eigenvalues 

We can use the asymptotic approximation obtained in QA.8P to reduce (25) to the simpler form 


C ° S ( 2 ^) (Jl (2 ^ } + ° (A)) + Sin ( 2 ^) (Fl + ° (A)) = °' 


This in turn can be rewritten as 


Applying Lemma [8] with 


k = 


ton (WL) = - Jl (Vyj + Q(A). 

V2 vx) Yi^Vr) 


277 


we obtain, after a minor effort, 


7 2 7 T 2 


An = ^^(nvr- 0 Go ( 7 )) + 0(n 4 ) 

2 


-*?( 1 + l ; eG » (7)+0< "" 2) ) 

= ir-( 1 -^ 0R ' h) + o{n ~ 2) ) 

as n —>■ +00, thus proving the first part of Theorem [4] 


(30) 


(31) 


(32) 

(33) 
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5.4 Approximation of Neumann eigenvalues 


The analysis for Neumann eigenvalues is slightly more complicated. Again we can use (A.8) to 

P( 7 , A) 


reduce (26) to 


tanp^ 

V 2 v/A 


Q( 7 , A) 


(34) 


where 


P(7, A) := - 8 )(\/A + 1) J\ (2^7) + (H^A - 87 X 2 VA - 7 ) J 2 ( 2 ^/ 7 ) 

- 8 V 7 ^A( 2 ^A - 7 )h (2^7) + 0(A 3/2 ), 

<3(7, A) := V 7 (sVA - 8 )(a/A + 1)11 (2^7) + (H^A - 87)(2v / A - 7 )F 2 (2^) 

- 8V7^A(2^A - 7 )y 3 (2^/7) + 0(A 3 / 2 ). 


Applying once again Lemma [8] with 

K = iTx’ Oo(7) = 


^(7,0) 

< 3 ( 7 , 0 ) 


Jo ( 2 y/ 7 ) 
Ao (2^7) 


-^0(7), 


we quickly arrive at 

Xn = h{ 1+ ^ 0Go(7)+ ° (n_2) ) = h ( x _ ^ 0%(7)+ ° (n_2) ) (35) 

as n —> 00, thus proving the second part of Theorem [4] 


6 Proof of the asymptotic results of the non-self-adjoint operator 


6.1 Eigenvalues and the Jost solutions 

Lemma 9. The eigenvalues of ([ 5 ]) are the zeroes of the determinant 

M{n) = = <^(7, 0) <£_ m (7, 0) + ¥>'_„(■ 7, 0) ^(7, 0). (36) 

Proof. Suppose that /i £ C is an eigenvalue of Ay, and that fj fl (x) G L 2 (M) is a corresponding 
eigenfunction. Then ^ solves the differential equation 

d 2 7 

~~ 1 + \ X \ 9 ^ X) = sign ( x )AS^( x )- 

If g± denote the restrictions of g M on M + and R_, then by integrability we must have 


g + (x) = C+tpui 7 ,x), 5 -(-x) = C'_< y o_ A1 (7,x), 


a; G 




where <^( 7 , x) is the Jost solution (24). 

As an eigenfunction should be continuously differentiable at zero, we obtain 


C+Mt, 0) - C'-<4’-m(7,0) = 0 , 

C+<4(7,0) + C-(p'_n(7, 0) = 0, 


which has a non-trivial solution if and only if M 7 (ju) = 0. 


□ 


Remark 10. (a) It follows from (24) that if g is real, either or is not square integrable, 

and therefore Ay cannot have real eigenvalues. 
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(b) By |BeTr, Proposition 4.6] one can instead look for the eigenvalues of ([3]) as the zeroes of the 
m-fu notion 

m („) = Ah.°> + 0) 

7 fctT.O) (P-/i(7.0) 

The use of half-line m-functions is natural and has been already suggested elsewhere, and 
described in great generality for indefinite Sturm-Liovuille problems with turning point at 0 in 
fKiTF] (see also references therein). 

(c) In what follows we assume that /j is in the upper half plane C + and look for the eigenvalues 
on the first quadrant. The final result will follow by symmetry (see Proposition 0 a ) and 
Proposition [ljf )). 



6.2 The determinant 


functions to rewrite (36) as 


We can use (24) and the known relations [ OlLoBoCl §13.3] between Kummer hypergeometric 

(' yV^M + 2ii)u( - —=,-l; 2 t^=, 0; 2^]i) 


M{li) = 


7 2 y/~ii e \/m 


8^ 5 /2 


2\/M 


+ 2// (y/—n + y/~ji — 2) u(^ 


2>/ W 


, 0 ; 2 % /=^)[/(- 2 ^=, 0 ; 2 ,/a) 


To find approximate solutions of M(a) = 0, we use the asymptotic formula (A.8). 
Let us define for brevity 


jv '■— Ji '( 2 ^/ 7 ), 
77r 


/T_l_ := cos 
£ 4 . := sin 
TV := tan 


Vv 

77r 


vW’ 

( 77 T \ 

IW’ 


:= cosh 

S7 := sinh 

T_ := tanh 


’ 

' 77T \ 

vW’ 

r 77T \ 

V 2 \/a/ 


(38) 


With these abbreviations and with the use of asymptotic formulae (A.8), equation (36) becomes 


7 


— r [-^- + 1 ] r . — 
Va \2^J1 ) \2 ^ 

Hi Va 


+ 1 


1 - 


) ( 2 a -wVa) + yi K ~) 


1 + 

+ i( 1 + 


87 

5Va 

87 


(j2 5+ +h/2 -S'-) + 


Va 

87 

x/aG^s 1 - - ij3<s , +) 


V7 


) C71-S+ + iyiS_) 


(2a - 7\/a) 


(39) 


(l _ g) (VA - (1 + i)) v/A (87 - 5^/a) (ji*4 + ViK-) 


Vi 


= 0 , 
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where we have dropped the lower order terms. 

Simplifying, writing S± = T±K±, and collecting terms in K±, we get 

K+{(1 + i)jl - (1 + i)) (64i7 2 - 40(1 + 1 ) 7 ^ + 25/x) 

+ T-yi (- 8 ry 2 + (5 + 16i)7^//I - 10/x) (j 2 (87 - lly/Ji) + 8 j 3 v / 7 ^u) 

- ji - 16ij 3 \/7v/M (4(1. + i)7 2 - 57\/m + 5(1 - i)/x) 

+ 2j 2 (64 7 3 - 128(1 - i) 7 2 ^ - 135i 7 /z + 55(1 + i)/x 3/2 ) 

+ T_ (87-5v// 7 ) ((1 + i)V 7 (v^~ (1 + i)) (87 + 5^) yi 
+ (7 + 2i v / ii)( 8 V7v/^2/3 + y2(-ll^ + 8i7))) } 

- K-[yi S'hy/ly/Ji' (- 8 i 7 2 + (16 + 51)7-^ - 10//) (40) 

+ 32 (647 s - (40 - 216i)7 2 ^/p - (176 + 135i)7/x + 110// 3/2 ) 

+ (1 + i)T_ {y/Ji - (1 + i)) (Wy 2 - 40(1 - 1 ) 7 ^ - 25ipt) yi 
+ 8(1 - i)-v/7\/W3 (4(1 + i)7 2 - 57^ + 5(1 - i)//) 

+ V 2 (64(1 + i) 7 3 - 2567 2 - v //z + 135(1 - 1 ) 7 // + 110i// 3 / 2 ) ) 

+ it (87 + 5iy^t) ((k - 2 V / / 1 ) (y 2 (87 - 11 y/y) + 8 ^/ 7 V//y 3 ) 

+ (1 - i)\/7(v 7 ^ -(! + !)) (87- 5 ^/ 1 ) i/i)} = 0. 

In what follows, we essentially replicate the reasoning in Lemma [8] but working to a higher 
accuracy. Introduce in the equation (40) the ansatz 




= A + v\ 3/2 + 7 A 2 , 


(41) 


where A is, as before, the absolute value of an eigenvalue of the self-adjoint operator with either 
Dirichlet or Neumann boundary conditions. Now replace back Ji_, T_ and K + with the corre¬ 
sponding expressions. The next step — expanding the left-hand side of the resulting equation in the 
Taylor series with respect to A around zero, — is delicate. 

First of all, observe that as A —> 0 


cos 


77T 


2^J \ + 1 /A 3 / 2 + i]\ 2 / 


and similarly 


( 77T 
= cos —= 
\2a/A 


H 




= COS 




/ 77r \ / 77T / ( 3 u z 

I \ „„„ I I v + ?? _ 


VA 


+ sin 


2 VA/ 1 " (t (- (- ¥)^ 


77T 


+ O(A), 


sin 


77T 


2\J A + z^A 3 / 2 + r/X 2 1 


= sin (^) cos (f 7\ 


77T A . / 77T / / 3z/ 2 \ 

- cosi 2 vt) sm Ur ("-t 


H) 


+ 0(A)- 
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We want to derive a similar expansion for tanh 


7 7r ] 

2 \J\+u\ 3 / 2 +r/\ 2 J 


. We use 


tanh(ti — t 2 ) 


sinh(f 2 ) + cosh(t 2 ) tanh(fi) 
cosh(f 2 ) — sinh(t 2 ) tanh(fi) 


with t\ := ^ and t 2 := Y (y + (rj - v/A + 0(A)). 

As tanh(7T7/2\/A) = 1 for A -7 0 modulo exponentially small terms, we get (again up to 
exponentially small errors) 


tanh ( — 77F - 

\2y A + 1 /A 3 / 2 + r]X 2 

and ( [40] ) reduces to an equation involving only tan ( 7 ^), tan [y + (y — ?,l f ) \/a)) and pow¬ 
ers of v/A. This is still, however, very hard to control. 


sinh(t2) + cosh(t2) 
cosh(t 2 ) — sinh^) 


6.3 Complex eigenvalue curves 

We can now use our knowledge of the self-adjoint problem (see (31) and fl34])) to simplify (40) 
further. For definiteness, suppose that —A is an eigenvalue of T D . 

Using the approximate identity (31), obtained for the eigenvalues of the Dirichlet self-adjoint 
problem on the half line, we can reduce the already simplified (40) to an equation involving only 


tan (y + (y — \/a)) an< ^ P owers °f '/A- Collecting the tangent terms, after some cu 


m- 


bersome but straightforward simplifications we arrive at 


where 



£(A,7) 
Q(A, 7 ) 


(42) 


P{ A, 7) = -(4 + 4i) 7 + (12 + 7i)VA, 

Q{ A, 7 ) = i( 4 + 4i)7 + (7 + 12i)v/A 

+ 2vr (4/y ((1 + 1)7 - Sv^) (joji + yoyi) + 7\/X (jf + Vi)) , 

and we have dropped terms of order 0(A). We have used the following standard relations [ OIL 0 B 0 CI 
Ch. 13] for the Bessel functions, 


Jn+i(2y/j)Yn(2y/j) - J n (2^y)Y n+ i(2^Y) = 

J i( 2 V^f) - ViM 2 \/i) = 

*1(27/7) - 7/7*2(27/7) = 7/7*0(27/7), 


in the simplifications. 

Note that right-hand side of (42) does not depend on u. We can now invert the tangent and 
solve with respect to is, to find the coefficient for the A 3 / 2 term in ©■ Expanding in A around 0 
and taking the leading term we get 


is =: T ( 7 ) = — arctan 

7T7 


1 


v i - 2V7 (</o(2V7) Jl{ 2 y/l) + 4o(2V7) *1(2^7)) ) ' 

The contribution of the other terms in the expansion then forms a part of ?y, which we drop. 


(43) 
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Remark 11. We can repeat the same procedure using the relation (34) for the absolute value of a 
Neumann eigenvalue A as a starting point. In this case we obtain 


v =: T + ( 7 ) = v = — arctan - 
7T7 V i 


1 


+ 27Ty / 7 (Jo(2V7) + Y 0 (2^y) n(2V7)) 7 ' 


(44) 


9(7) 


Remark 12. One can use standard asymptotic formulas for Bessel functions (see [ OIL 0 B 0 CI , Chapter 
10 ]) to observe 

_l 2 iM if 7^0 

J_ _ 3cos ( 4 ^) 

Moreover, again using standard estimates and properties of Bessel functions and their zeroes one 
can observe that 

| J 0 (x)|, | |Y o 0c)|, l^tOzOI < Wf for x - 2 ’ 

and the four functions are monotone for x < 2 (Jo increasing and bounded by 1 , the other ones 
decreasing, in particular the Y are unbounded). In particular, for any 7 e (0,1) neither T~( 7 ) and 
T + ( 7 ) nor their real and imaginary parts, vanish. Moreover 7 ^( 7 ) define curves in the complex 
plane that diverge as 7 —»• 0 and converge to 0 as 7 —»• +00. 

With q{ 7 ) as in © we get fll3] ). Then the identity ImT ( 7 ) = ImT + ( 7 ) follows immediately 
using the standard relation between arctan and log (see Remark [7](a)), thus proving part (i) of 
Theorem [ 6 ] 

Proving part (ii) of Theorem [ 6 ] requires some extra work. First of all observe that, up to the 
errors of order 0 (n~ 4 ), we have 


X n~ X n=— 2 ^(® R ,h)-®Roh)) 


7 


47rn 3 
,,2 


7 


47rn 3 


ilog 

ilog 


’ (Ji(2/y) - ili(2^7))( J 0 (2^) + iY 0 (2^)) \ ' 

(J 1 (2V7) + iF 1 (2^7))(Jo(2V7)-i>o(2V7))| 

g( 7) +i N 

9 ( 7 ) - i 


where Rki'j) are defined in Theorem [ 4 ] and the first identity follows from the fact that zeroes of Yq 
and Y\ are interlaced. Similarly 


r-(7)-T77) = ^ log (|7^|) 


29 ( 7 ) + 2i/ ' 

Despite the appearance of the complex unity i in the above formulae, all these expressions are in 
fact real! 

To show that the two asymptotic formulae (12) for fx n coincide up to the lower order terms, we 
use Theorem [4] to obtain, as n -7 00, 


(A? + T -( 7 )(\°) 3/2 ) - (A? + T+( 7 )(A^) 3 / 2 ) + O 


rr 


7 2 i 

47rn 3 

7 2 i 

47 TO 3 


log 


log 


g( 7) + i 
9(7) - i 
g( 7) + i 
g(7) - i 


2i , 

H-log 

7T7 


2g(7) + 2i 


- log 


2 ( 7 ( 7 ) - 21 

g( 7) +1 
g(7) - 1 


-) f—) 

i ) l 4n 2 / 


3/2 


+Oh 


rr 


+ o[- i 


= o 


n 


4 / ’ 


thus concluding the proof. 
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7 Generalizations and other remarks 



Figure 5: Approximated eigenvalues of y 4 ( 7 + , 7 _) for 7 _ = 1.5, 7 + = 5. 

The procedure used to prove Theorem [6] can be repeated in a completely similar way to obtain 
a result for the operator 


A( 7+ ,7_) = JT V , V{x) 


7+ 

i+l*l 

7- 

1 + 1*1 


if x > 0 
if x < 0 ’ 


7+! 7- e M+. 


In this case the m-function is of the form 


M( A) 


<7 + (°) ! ^,7-(°) - 

¥’/i, 7 +( 0 ) 


The curves in the upper (resp. lower) half plane are no more symmetric w.r.t. iM, however for the 
left quadrants and right quadrants we can extend Theorem [6] The only difference is that now the 
Y” and Y + are now functions of both 7 + and 7 _. 

Let z/, 77 E M+. Set 


f-M ■= + r„(2^)F l( 2^))) 

- iM 2 Vv)Ji( 2 Vv) + Y 0 (2y/fj)Y 1 (2y/fj)) , 

f+M ■= X ( ! + ^ + y.(V?)y.(Vg))) 

+ TTy/rj (J 0 (2 v Yj)J 1 (2^/rj) + Y 0 (2y/rj)Yi(2y/rj)). 


Then the two factor multiplying the term Re/ 7 3 / 2 are given by 


T (7+ >7-) 
T + (7+>7-) 


£ f 7- 1 arctan(l//_( 7 +, 7 _)) 
7T [ 7+ 1 arctan(l//_( 7 _, 7 +)) 

£ f 7 Z 1 arctan(l//+( 7 + , 7 _)) 
vr [ 7- 1 arctan(l//+( 7 _, 7 + )) 


if Re n > 0 
if Re n < 0 

if Re n > 0 
if Re n < 0 


One can immediately see that the asymmetry appearing w.r.t. iM is reflected in the asymmetric 
dependence on 7 + and 7 _. 
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It is interesting to observe that for Re/i > 0 the effect of 7 _ is much stronger than the one 
of 7 + (the latter appears only in the cotangent term, its contribution is bounded, while the former 
additionally appears as an inverse prefactor). The situation is opposite when Re fi < 0. 

The expressions for Y T ( 7 _, 7 + ) are more involved than the ones for Y T ( 7 ) but, as expected, 
they simplify to (43) and (44) for 7 + = 7 _. As that case, it is possible to use the standard results 
on Bessel functions to show that the two constants have non-zero real and imaginary part for any 
7± > 0. 


To answer the general question posed in Bel3 j for a wider class of potentials one would need 
good estimates of the Jost functions in a complex half ball containing the origin and the positive 
and negative real axis. To our knowledge, the best result of this kind is contained in a paper by 
Yafaev [Yaj . In that work, however, the author needed to exclude two cones containing the real axis 
for his estimates to hold. Additionally he could get only the first term in the asymptotic expansion, 
whereas for our result we would need at least the first two. 


A Uniform asymptotic expansion of Kummer Hypergeometric func¬ 
tions 

We need to approximate 

u ( ~ ^=,c; 2 % / = a), C/(-^=,c; 2^), cG{0,-l}, (A.l) 

in the limit /i —>• 0. We use the theory developed in [Te. Chapter 27]. By formula |Tel (27.4.85)], 
as a —> 00 , and with az bounded and Re(az) > 0, 

/ 00 4 00 D \ 

U(-a,c ; az) ~ /5 1 ~ c T(a+ l)e^ C c ^(Q £ + /?C c _ 2 (C) £ , ( A -2) 

V n =0 a n=0 a ) 

where A n and B n are defined by an iterative procedure, ( = 2 (3a and 

C„(Q = cos(7ra) J V (Q + sin(7ra)lY(C)- (A.3) 

On can immediately see that in our case 

a = a± = , z = z± = - T . (A.4) 

2 v /± 7 t a z ± 

We will drop the subscript ± for the rest of the discussion. 

Additionally we need an expression for j3. This is defined in [Tel (27.4.36)] as 

(3 = + sinh(u? 0 )) 

where wq = 2 arcsinh(^ sfz), see [Te, (27.4.33)]. Equation [Te, (27.4.52)] gives an asymptotic 
expansion for /3 as z —> 0 : 

P 2 = z+^z 2 + 0(z 3 ). 

By Taylor expansion we get 

P = \fz\j 1 + Y 2 Z + °( z2 ) = ( x + + °(- 2 )) • ( A - 5 ) 

Therefore 

C = 2 ^a = 2^/7 ^1 + " + 0 (a 4 )^j 
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Figure 6: Plot of real part (left) and imaginary part (right) of U (black) and 

its approximation given by QA.8P (dashed red) for small values of /i and 7 = 2.5. 


and Re(ac:) = Re( 7 /a). 

Observe that for 7 € M+, Re(ac:) > 0 iff Re(a) > 0. 

The coefficients Aq and Bq also have explicit expressions that can be derived using some sym¬ 
metry properties and L’Hopital rule, see [Tel (27.4.74)]): 


An = 


B 0 = 


( p \ 

| \ 

V2 sin(0 )) 

( P \ 

1 V 

\2sin(0) J 

1 V 


— tan 6 cos (c6) 


2 sin(d?) 

~V 


where 9 = — ^i^o- 

The computation of A n and B n for n > 0 is quite involved, however we will need only A\. One 
can exploit the procedure to compute Aq and Bq, and the recursive definition of the coefficients to 
get a Taylor approximation in negative powers of a for c E {0, —1}. We get 


if c = 0, TLq = 1 + 0(a~ 2 ), A? = -- + 0(a- 2 ), (3B° 0 = 0, 

Id 


(A.6) 


if c= -1, Tg 1 = 1 + 0(a -4 ), A^ = -^ + 0(a~ 2 ), pBp 1 = + 0(a~ 3 ). (A.7) 


With these, (A.2) can be re-written 


U (~a,c; 7 

V a 

where 


1 — c 

y/l^ 2 


T(a+l)e^ (c c -i(a, 7 )(^ + ^) + C' c - 2 (a, 7 )^ 0 C + O(a- 2 )) (A. 8 ) 


C u (a, 7 ) := cos( 7 ra)J„( 2 x / 7 ) + sin( 7 ra)y],( 2 v / 7 ) 


(A.9) 


and Aq, A\ and Bq are obtained dropping the error term in the appropriate coefficient in (A.6) and 

63- 

Remark A.l. Here the error is in fact O^/a 2 ), we may thus expect the improvement in the precision 
of the asymptotics when 7 > 1. 

Remark A.2 (Validity of the expansion). If we define 

h = (3 + vri + yJ(P + 7ri) 2 - (3 2 , 


then the asymptotic formula (A.2) is valid for 


7T 7T 

- arg t 1 — - + 0 < arg a < arg h + — — 0 , 


(A.10) 


and the same applies to z (see [Te] Chapter 26.4.2]). 
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In our case 


ti ~ 27ri + 2/3 + 0(/3 2 ). 

For \a\ 1, argti is in the upper complex half plane. In particular this allows a and z to be in the 

closure of the first and fourth quadrant. 
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